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EXECUTIVE  SUMMARY 


We  develop  the  fundamental  exact  analytical  and  computational  model 
required  to  study  the  transmission  of  incident  plane  sound  waves  into 
submerged  elastic  cylindrical  shells  subjected  to  arbitrary  loads  on  their 
outer  surface.  We  use  the  superposition  principle  in  this  linear  problem  to 
separate  the  contributions  to  the  internally  transmitted  field  caused  by  the 
incident  wave  from  that  of  the  surface  load.  This  basic  model  uses  the  exact 
(2-D)  formulation  of  elastodynamics  to  describe  the  shell  motions,  and  that  of 
general  linear  acoustics  to  describe  the  wave  motion  in  the  inner  and  outer 
fluids.  We  display  the  isobaric  contours  of  the  internally  transmitted 
pressure  fields,  exhibiting  their  caustics  and  their  progressive  development 
as  the  frequency  is  increased  within  the  band  0  <  k3a  <  10.  The  contour  plots 
are  generated  for  all  the  possible  combinations  of  loading  and  insonifi cation 
conditions  with  a  view  toward  computing  an  advantage  ratio  PR  which  measures 
the  gains  that  would  result  from  sensing  the  internally  transmitted  field, 
rather  than  the  -field  external  to  the  shell.  This  advantage  ratio  is  shown 
always  to  be  greater  than  unity,  and  in  many  cases,  much  greater  than  unity. 
The  effects  of  two  basically  different  types  of  loads  are  analyzed  via  the 
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frequencies.  We  observe  an  encouraging  overall  permanence  of  the  caustic 
locations,  at  fixed  frequencies,  as  the  types  of  surface  loads  are  varied. 

The  present  investigation  of  sound  transmission  into  arbitrarily  loaded  shells 
clearly  demonstrates  the  filtering  behavior  of  the  shell  in  the  frequency 
domain,  and  its  focusing  action  in  space.  Hence,  for  moderate  surface  loads, 
the  air-filled  submerged  shell  analyzed  here  acts  as  a  very  effective  sound 
concentrator,  particularly  near  the  resonances  of  its  monopole  mode. 

We  have  carried  on  a  detailed  computational  study  in  the  final  sections 
of  this  report  to  quantitatively  determine  the  effect  of  shell  thickness, 
stiffness,  and  general  material  composition,  on  the  sound  fields  transmitted 
into  the  loaded  shells.  This  analysis  makes  the  work  valid  for  various  metals 
and  for  shell  relati ve-thickness  between  1  percent  and  20  percent,  which  is  a 
range  covering  most  practical  cases.  This  final  section  shows  that  for  all 
thicknesses  and  compositions,  the  large  values  of  the  advantage  ratio  PR, 
still  show  large  gains  in  sensing  internal  rather  than  external  fields  in  all 
cases  considered,  provided  the  forcing  (flow)  function  is  below  acceptable 
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FOREWORD 

This  report  presents  a  detailed  analytical  and  computational  study  of  the 
sound  transmission  process  that  takes  place  when  an  elastic  cylindrical  shell 
submerged  in  water  is  subjected  to  the  double  excitation  caused  by  an  incident 
acoustic  wave  and  by  the  action  of  surface  forcing  functions,  of  various 
essentially  different  types,  applied  along  the  periphery  of  the  shell.  We 
constructed  a  novel  analytic  methodology  to  evaluate  the  sound  pressure  fields 
inside  the  shell  caused  by  all  possible  combinations  of  insonifi cation  and 
surface  excitation  conditions.  This  methodology  is  used  to  compute  an 
advantage  ratio  PR  which  measures  the  gains  that  would  result  from  sensing 
internally-transmitted  fields,  rather  than  fields  external  to  the  shell. 

This  work  was  partly  supported  by  Code  432  of  the  Office  of  Naval 
Research  and  partly  by  the  Independent  Research  Board  of  the  Naval  Surface 
Weapons  Center  during  FY83  and  FY84. 
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THE  CYLINDRICAL  SHELL  WITH  LOAD  #2  ACTING  ON  PORTION  OF  ITS 
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INTERNALLY  TRANSMITTED  ISOBARIC  CONTOUR  PLOTS  FOR  AN  AIR-FILLED 
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INTERNALLY  TRANSMITTED  ISOBARIC  CONTOUR  PLOTS  FOR  A  SUBMERGED, 
AIR-FILLED,  CYLINDRICAL  SHELL,  LOADED  WITH  LOAD  #1  (UPPER  PLOT), 
OR  WITH  LOAD  #2  (LOWER  PLOT),  BOTH  FOR  k3a  =  5.  NO  INCIDENT 

WAVE  INTERACTS  WITH  THE  SHELL  HERE  . 
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INTERIOR  PRESSURE  LEVELS,  |p3/Pq|,  VERSUS  k3a,  AT  FIVE  (OF  ELEV¬ 
EN)  LOCATIONS  INSIDE  AN  UNLOADED  CYLINDRICAL  SHELL  IN  WATER 
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INTERIOR  PRESSURE  LEVELS,  1 P3/P0 1 .  VERSUS  k3a,  AT  FIVE  (OF  ELEV¬ 
EN)  POINTS  INSIDE  A  SUBMERGED  CYLINDRICAL  SHELL  LOADED  WITH  LOAD 
#1 ,  AND  INSONIFIED  BY  AN  INCIDENT  PLANE  WAVE  . 
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SECTION  1 
INTRODUCTION 


The  scattering  of  sound  waves  incident  on  cylindrical  elastic  shells  has 
traditionally  been  studied  in  the  outer  fluid  surrounding  the  shells, 
particularly  at  far  away  distances.  This  is  a  consequence  of  the  fact  that 
this  is  the  operating  region  of  most  sensing  devices.  The  literature  on  this 
type  of  sound/structure  interaction  problems  is  quite  voluminous1"12  and  it  is 
not  possible  to  cover  it  extensively  here.  Some  authors  have  treated  a  single 
bare  shell. 1,7,9  Others  have  studied  layered  concentric  shells,5,6,0  either 
all  elastic11  or  some  elastic  and  some  viscoelastic.6,10  Some  have  described 
the  shell (s)  motions  by  the  exact  equations  of  linear  two-dimensional 
elastodynamics  (viscoelastodynamics)  while  others  have  described  the  shell 
responses  by  means  of  the  equations  of  some  approximate  shell  theory.2,3  Some 
have  considered  the  same  fluid  inside  and  outside  the  shell,  while  others 
have  studied  the  shell  separating  dissimilar  fluids.5,12  Sometimes  the 
analysis  has  focused  on  techniques  for  the  low-frequency  response7  while  other 
authors  have  emphasized  methods  particularly  useful  at  high  frequencies.1,6 
Some  studies  have  shown  connections  between  the  scattering  and  the  radiation 
of  sound.  In  most  cases1"12  the  incident  waves  are  considered  plane  as  they 
would  be  if  they  emerged  from  a  distant  transducer,  although  some  isolated 
cases  of  spherical  wave  incidences  have  been  treated.13"15  For  incident  plane 
waves,  almost  all  the  literature  has  dealt  with  the  scattered  waveforms  in 
the  outer  medium,  far  away  from  the  shell.  We  know  of  hardly  any  studies  of 
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backscattered  near-fields.  In  all  cases  that  we  have  ever  encountered,  the 


shells  are  usually  immersed  in  fluids,  but  are  otherwise  unloaded. 


In  the  present  paper  we  will  study  the  transmission  of  sound  energy  into 


the  fluid  medium  inside  a  single,  elastic,  and  loaded  cylindrical  shell 


submerged  in  a  dissimilar  fluid.  We  are  not  aware  of  any  published  work  on 


these  aspects  of  the  interaction.  The  incident  acoustic  wave  is  plane  and  it 


sets  the  shell  into  vibration  as  it  falls  upon  it.  The  shell  is  free  to 


translate  by  radiation  pressure  as  the  wave  hits  it.  The  shell  is  of  a  great 


length,  and  the  incidence  is  at  broadside  aspects.  In  this  common  two¬ 


dimensional  situation,  the  shell  motions  are  governed  by  the  corresponding  z- 


independent  equations  of  linear  elastodynamics.  The  shell  is  not  only  excited 


by  the  incident  plane  harmonic  wave,  but  it  is  also  subjected  to  an  arbitrary 


loading  on  its  outer  surface  caused  perhaps  (although  the  cause  is  not 


important  here)  by  a  relative  motion  of  the  outer  fluid  past  the  shell.  So, 


the  shell  is  not  only  hydrostatically  fluid-loaded  by  the  simple  fact  that  it 


is  submerged  in  a  fluid,  but  also  because  there  is  an  externally-generated 


arbitrary  load  acting  on  its  outer  surface.  We  wish  to  systematically  study 


the  internal  sound  fields  developed  inside  the  shell  as  a  result  of  the  double 


action  of  the  incident  plane  wave  of  amplitude  Pq  and  of  the  arbitrary  surface 
load  of  amplitude  p£.  It  becomes  immediately  obvious  that  the  innumerable 


modal  shell  resonances  that  are  excited  by  this  double  action,  will  induce 


various  caustics  in  the  internally  transmitted  sound  pressure  fields,  and  that 


at  these  selected  points,  the  internal  pressures  could  be  very  pronounced.  We 


will  see  how  these  intuitive  empiricisms  actually  emerge  and  are  substantiated 


by  the  following  model  and  its  analytic  predictions. 


raw 


w 


5X* 


NSWC  TR  85-42 


SECTION  2 

THEORETICAL  BACKGROUND 


The  theory  has  already  been  developed  for  bare  and  coated  but  unloaded 
cylindrical  shells. 1,6  We  briefly  outline  the  basics®  for  completeness. 
Navier  equations  for  the  (infinite)  cylindrical  geometry  are 


vre  \ 


“r  2  ,u8  >  . 

7-7-35- +  — i 

r  r 


3A  .  2 
3r  ks  ur 


v2au  Zt  t  iJ_ii  i  |i  +  k  2„  o 

re  9  7  7  ,9  y  r  ae  s  e 


where  the  Laplacian  and  the  dilatation  are,  respectively. 


„  2  2 
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7  re  ■  — 2  *  7  a?  — s 
r0  3r  r  3r  r  30 


A  .  ^r  .  1  %  +  V 
A  "  3r  r  30  r  » 


2  2  2  2 

where  k$  =  po>  /u,  k^  =  pw  /(A  +  2y),  and  where  time-dependence  of  the  form 
exp(-iwt)  has  been  assumed  throughout. 

We  have  shown  earlier16  that  Equations  (1)  are  automatically  satisfied  if 
two  scalar  potentials  are  introduced  which  satisfy  the  two  scalar  Helmholtz 


equations 


<'r 0  +  k<|2>*  "  °> 


<’r!  *  ks\  ’  0  ’ 
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All  the  relevant  displacement  and  stress  components  can  be  obtained  from 


and  by  means  of  the  following  formulas:^5 


ur  3r  r  36 


u 

ue  r  36  3r 


2  a2.  a  i 
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3r36  v r 
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Solving  the  uncoupled  Equations  (3)  subject  to  suitable  boundary  conditions, 
which  are  set-up  by  means  of  Equations  (4)  and  (5),  one  is  really  solving  the 
coupled  Equations  (1).  Once  the  fully  determined  4>  and  t|»z  are  found,  then  the 
same  Equations  (4)  and  (5)  yield  the  displacement  and  stress  fields 
everywhere. 

NORMAL  MODE  SOLUTION  FOR  THE  UNLOADED  SHELL 

This  solution  can  be  obtained,  for  an  unloaded  shell,  from  earlier 
published  work.5  In  the  three  media  one  has, 

Pi(r,e,t)  =  pQe  “iwt  l  in  en  rJn(kir)  +  bnHn(1)(kir)]  cosne  (6) 


for  the  outer  fluid  (medium  1),  together  with 


(1)_  J_ 

ur  2  3r  * 

Pl<o 


In  the  shell  material  (medium  2),  the  elastic  potentials  are 


♦2  -  P0e'1Bt  j/  en  C<>n  J„<kd2r)  *  hn  Vkd2r)]  cos"9 


*z2  ’  "o5’01  <kS/>  *  Vn(V>5  s,"ne> 


(8) 


where  k.  =  to/c.  and  k  =  w/c  .  Finally,  the  pressure  in  the  inner  fluid 

02  d2  S2  S2 

(medium  3)  is 


P3(r,e,t)  •  Poe“1a>t  l  inenqnJn  (k3r)  cosn6> 

0  n=0  n  n  n 


(9) 


where  k3  *  U1/C3  and  the  six  sets  of  constants  bn,  gn,  hn,  tp,  mn,  and  qn  are 
determined  from  the  boundary  conditions,  which  for  the  unloaded  shell  are. 
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rr 
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re 


(2) 


(10) 


and  where  the  Newmann  symbol  is  en  =  2  -  6^  .  We  will  treat  the  case  of  the 
loaded  shell  later.  Substituting  Equations  (6) -( 9)  into  (10)  yields  the 
following  vector/matrix  linear  relation 


°v  V  • 


(id 


where  xn  =  (b^,  gn,  hn,  *n,  mn,  qn)T  and  X*  =  (A^,  A2*,  0,  0,  0,  0)T,  and 
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D  Is  a  6  x  6  matrix  that  can  be  easily  obtained  from  the  10  x  10  matrix  found 
earlier^  for  coated  shells.  The  result  is. 


dll 

d12 

d13 

d14 

d15 

0 

d21 

d22 

d23 

d24 

d25 

0 

Dn  = 

0 

d32 

d33 

d34 

d35 

0 

(12) 

0 

d86 

d87 

d88 

d89 

d8,10 

0 

d96 

d97 

d98 

d99 

d9,  10 

0 

d10 ,6 

d10. 

7  d10,8 

d10,9 

0 

where  these  elements  are 

C  O 

all  listed  before.  *  We  could  relabel  them  as 

fol lows : 

dll 

d12 

d13 

d14  d15 

0 

d21 

d22 

d23 

d24  d25 

0 

O 

13 

II 

0 

d32 

d33 

d34  d35 

0 

(12') 

0 

d42 

d43 

d44  d45 

d46 

0 

d52 

d53 

d54  d55 

d56 

0 

d62 

d63 

d64  d65 

0 

so  that  the  indices  run  from  1  to  5,  as  they  should  in  this  (uncoated)  case. 
Equation  (11)  can  be  solved  for  any  of  the  six  sets  of  constants  it  contains, 
by  Cramer's  rule.  For  example,  the  set  of  constants  qn  is  given  by  the 
quotient 
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where  On  is  as  above,  and  Qn  is  like  Dn  but  with  its  sixth  column  replaced  by 


the  column  vector  X 


%  '  <’/Dn) 


i  .e. 

11 

> 

d12 

d13 

d14 

d15 

Aj* 

21 

d22 

d23 

d24 

d25 

i 

d32 

d33 

d34 

d35 

0 

i 

d42 

d43 

d44 

d45 
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i 

d52 

d53 

d54 

d55 

0 

i 

d62 

d63 

d64 

d65 

0 

which  completes  the  solution.  The  30  different  elements  contained  in  the 
numerator  and  denominator  of  this  equation  are  listed  in  Appendix  A  to  make 
this  work  self-contained,  although  they  could  be  obtained  with  relative  ease 
from  our  earlier  work.*’*®  Finally,  it  follows  from  Equation  (9)  that 


p"  -  ■  l  1  S  %  Jn  ,k3r> 

u  n=o 


cosne 


(14’) 


THE  LOADED  SHELL  CASE 

If  the  outer  surface  (i.e.,  r  =  a)  of  the  shell  is  additionally  loaded 
with  a  spatially  varying  (i.e.,  8  -  dependent)  excitation  pL(a,e,t),  then  the 
normal-mode  solution  given  above  continues  to  hold,  but  now  the  first  of  the 
boundary  conditions  changes.  The  first  boundary  condition  in  this  loaded  case 
is 

Trp2\=a  =  -  CPi(r,e,t)  +  pL(r,e,t)]r=a  ,  (15) 


where  p^  is  still  given  by  Equation  (6)  and  P[_  is  the  new,  arbitrary  surface 
load.  Let  this  load  be  harmonic  in  time  and  spatially  symmetric  about  the 
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direction  6  =  0  of  the  incident  wave.  The  load  is  then  an  even  function,  and 
by  means  of  a  Fourier  (cosine)  transform  we  can  write 

p,(a,e,t)  =  p2e“lu)t  l  ine  cosne  [— \  p,  (e)cosnede],  (16) 

L  n  7T  ft  L 


where  p2  (*  PQ)  is  the  amplitude  of  this  load.  We  recall  that  the  incident 
plane  wave  is  given  by  the  expansion, ^ 


•  00 

Pinc(r,e,t)  =  p0exp  C i ( it ! «r  -  at)]  =  p0e'1a)t  l  i^J^lqr)  cosne.  (17) 

n=0 


In  view  of  the  fact  that  there  is  a  new  boundary  condition  at  r  =  a,  all 

the  sets  of  constants  (i.e.,b,c,...  q  )  for  this  loaded  case  will  be 

n  n  n 

different  from  before,  and  that  is  why  we  have  labeled  them  with  tildes. 

If  we  substitute  the  solutions  into  the  boundary  conditions  to  determine 
the  new  sets  of  constants,  and  we  concentrate  on  the  set  qn,  it  turns  out 
that  we  find  an  expression  very  similar  to  that  found  for  the  unloaded  case 
and  that  is  why  we  introduced  it  earlier.  The  result  now  is 
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where  Dn  is  still  given  in  Equation  ( 12 * )  and  the  above  determinant  is  almost 
identical  to  that  in  Equation  (14),  except  for  the  element  in  its  upper  right 
corner  which  in  Equation  (14)  was  simply,  (c/o  Appendix  A) 

A*=  -  (pi/p2)k$*  a2  Jn(M)  (19) 

while  now  in  Equation  (18),  with  the  tilde,  it  takes  the  form 

%i*  ■  Ai*  +  I  -  -  (pi/p2)k  2  a  [J  (kxa)  +  ~  ^  f  Pi (®)  cosnede],  (20) 
n  s2  n  Po  n  •'0  L 

Clearly,  the  portion  labeled  In,  containing  the  integral,  is  the  portion  that 
accounts  for  the  surface  loading  pl,  which  is  now  assumed  of  unit  amplitude 
(see  Figure  1)  since  its  amplitude  P2,  normalized  to  Pq  has  been  taken  out  of 
the  integral  and  placed  in  front  of  it  as  a  factor. 

For  the  internal  sound  pressure  levels  (SPL)  inside  the  shell  (i.e., 
medium  3),  Equation  (9)  can  be  rewritten  in  normalized  form  as  follows 

I  v"  qnlVlt3r)cosn0|  »  (21) 

where  qn,  for  the  loaded  shell,  is  given  in  Equation  (18)  with  Aj*  in  Equation 
(20).  This  is  the  key  result  for  the  systematic  study  of  the  internal  field 
of  the  loaded  shell  insonified  by  the  incident  plane  wave.  No  shell  theories 
have  been  introduced.  The  results  are  all  exact  under  the  conditions  stated. 
Note  that  qp  is  the  only  one  of  the  new  sets  of  constants  that  enters  the 
series  expression  for  the  (normalized)  pressure  field  inside  the  arbitrarily 
loaded  shell.  Most  of  the  following  sections  will  study  various  properties  of 
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the  solution  given  by  Equations  (21)  and  (20).  Before  the  advent  of  the 
present  computer  age,  it  would  have  been  unthinkable  to  evaluate  formulas  such 
as  these.  Note  that  coefficients  qp  and  qn,  from  Equations  (18)  and  (14), 
for  the  loaded  and  unloaded  cases,  respectively,  contain  all  the  necessary 
dependence  on  geometry,  on  material  composition  of  the  shell  and  of  its 
surrounding  fluids,  and  on  the  amplitudes  of  the  incident  wave  and  the  surface 
load,  as  well  as  on  the  spatial  distribution  of  the  surface  load.  All  these 
coefficients  also  depend  on  frequency  and  this  general  solution  is  valid  for 
the  entire  frequency  spectrum  including  its  low,  middle,  and  high  regions.  In 
view  of  the  usual  slow  convergence  of  series  solutions  containing  Bessel 
functions  such  as  that  in  Equation  (21),  it  may  be  computationally  expensive 
and/or  analytically  unfeasible  to  carry  on  meaningful  calculations  beyond  a 
value  of  k3a  of  about  15  or  20.  Beyond  this  value,  the  series  must  be 
conditioned  via  the  Watson-Sommerfeld  method^  or  other  similar  techniques,  not 
only  to  accelerate  the  series  convergence,  but  to  bring  out  the  physics.  At 
these  high  frequencies,  the  motivation  for  the  present  work  on  internal  fields 
within  shells  loses  its  urgency  and  it  will  not  be  pursued  here. 
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SECTION  3 

SUPERPOSITION  PROPERTIES  OF  THE  LOADED  SHELL  SOLUTION 

In  view  of  the  linearity  of  this  problem  we  can  study  the  effect  of  the 
incident  monochromatic  wave  on  the  internally  transmitted  field  separate  from 
the  effect  of  the  surface  load.  From  Equations  (18)  and  (20)  it  follows  that 

%  ■  %  +  %  <22> 


where  is  the  coefficient  associated  with  the  unloaded  shell  excited  only  by 
the  incident  plane  wave  already  given  in  Equation  (14),  and  q^  is  a  new 


coefficient  given  by  the  determinant  ratio 
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where 


Pi  2  2  P2  if* 

I  = - k  a  (— )  —  I  p,  (6)  cosnedQ. 

n  p2  s2  Po  *  J  o  L 


.  ,  -  _  .  m  V  1  •  •  •  •  •  »  •  «  »  •  ,  •  K  *.  *  '  V  •  *  *  1  '  »  *  ‘  ■  k  *  V*  »  ' 

•.V.‘  \VvV.\V  •  *  ••  •  •  •*  ■  '  O  v  ■.  ■w-v"  «.*  -  '  •  *  ■ v' 
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This  new  nondimensional  coefficient  represents  the  contribution  to  the 
internally  transmitted  field  caused  by  the  (isolated)  surface  load  in  the 
absence  of  the  incident  wave.  This  contribution  clearly  vanishes  for  unloaded 
shells.  The  determinant  Dn  is  still  as  given  in  Equation  ( 12 * ) . 

For  ease  of  calculation,  and  in  order  to  best  be  able  to  interprete  the 
steps  of  the  analysis  that  will  follow  in  this  section,  we  will  rewrite 
Equation  (22)  in  the  form 

P2  - 

q  =  q  +  ( — )  q  (25) 

Mn  Hn  'p0'  Mn  '  1 


where  qp  now  is  given  by 
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and  where 


-  PI  2  2  if 

I  = - k  a  —  |  p,  (e)  cosnede. 

n  P2  S2  *Ja  L 


We  have  simply  explicitly  pulled  out  the  factor  (Pz/pg)  from  Equation  (24) 
It  is  now  possible  to  return  to  Equation  (21)  which  gives  the  (normalized) 
internal  sound-pressure  level  of  the  loaded  shell,  and  analyze  its  two 
components  given  by  Equation  (25)  separately.  We  define 


V  h  s  A  v  %  (k3r>  cosn9  • 


>  A  A 
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as  the  contribution  to  the  SPL  caused  solely  by  the  incident  wave,  and 

P2  - 


P3 

P  =  |  — - 
Po 


l  Jn(k3r)  cosne 


n=0 


n  po  ^n  n 


(29) 


as  the  contribution  to  the  internal  SPL  solely  due  to  the  surface  load. 
Multiplication  through  by  (pg/P2)  yields. 


po  * 
P2  P 


A 

P3 

pT 

A 

I  O  qnJn(k3'*)cosne 
n=0  n  "  n 


(30) 


We  now  consider  the  ratio  of  the  contributions  of  Equations  (28)  and  (30)  to 
the  total  internal  SPL.  This  would  be  the  ratio  of  the  (duly  normalized) 
contribution  to  the  internally  transmitted  pressure  field  caused  by  the 
incident  wave  alone,  to  that  of  the  contribution  caused  by  the  surface  load 
alone.  This  resulting  ratio  PR  is. 


|  p  3/ P  0 

A 

A 

P3/P3 

nq„J „(  k 3r)  cosne 
n=0  n  ^n  n'  3 

A 

1  A 

P3/P2 

P0/P2 

CO  2 

M  A 

I  O  qnJn(k3r)  cosne 
n=0  n  n  n 

(31) 


This  interesting  result  shows  the  ratio  of  the  incident  signal  to  surface-load 
amplitudes,  when  they  are  computed  inside  the  shell  (numerator)  and  outside 

A 

the  shell  (denominator).  Coefficients  qn  ,  qp  are  given  in  Equations  (14)  and 
(26)  respectively.  If  the  surface  load  is  to  be  viewed  as  some  undesirable 
"noise"  that  corrupts  the  incident  waveform,  then  the  ratio  PR  can  be 
interpreted  as  the  "signal -to-noise  ratio"  sensed  after  transmission  into  the 
shell,  divided  by  the  signal -to-noise  ratio,  sensed  outside  of  it.  If  PR  >  1 
it  then  follows  that  there  would  be  an  advantage  in  sensing  these  fields 
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internally  rather  than  externally. 

Finally,  if  we  want  to  study  the  total  internal  pressure  field  containing 
the  contributions  from  both  incident  signal  and  the  arbitrary  surface  load, 
the  result  is  that  given  in  Equation  (21),  with  qn  given  by  Equation  (25). 
This  normalized  result  is, 

Jo  e"in|:q"  +  W  q"] J" (k3r)  cos"9  ’  (32) 

where  qn  and  qn  are  respectively  given  in  Equations  (14)  and  (26).  The 
result  computed  from  this  formula  would  not  differ  much  from  the  analogous 
result  for  an  unloaded  shell  given  in  Equation  (28)  if  the  effect  of  the 
surface  load  were  negligible.  This  obvious  criterion  will  serve  to  determine 
if  a  given  surface  load  is  negligible  or  not.  We  will  see  that  there  are  some 
types  of  surface  loads  that  affect  just  a  couple  of  the  shell's  normal  modes 
of  vibration  and,  thus,  have  a  relatively  small  disturbing  effect  on  the 
pressure  field  transmitted  into  the  shell  as  a  result  of  the  (dominant) 
incident-wave  contribution.  Other  types  of  surface  loads  affect  many,  if  not 
all,  of  the  shell's  normal  modes  of  vibration  and,  therefore,  have  a 
pronounced  distorting  effect  on  the  internal  fields  caused  by  the  incident 
wave  contribution. 

Before  showing  numerical  results  for  specific  loads  and  geometries,  we 
can  make  some  general  remarks  based  on  physical  intuition  and  common 
scientific  sense.  As  the  outer  shell  to  sound  wavelength-ratio  (a/X)  is 
increased,  a  more  and  more  complex  sound  field  will  be  transmitted  into  the 
shell's  interior  and  families  of  caustics  will  inevitably  be  developed 
internally.  If  we  plot  the  internal  isobars  (i.e.,  the  loci  of  points 
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exhibiting  the  same  value  of  the  internal  sound  pressure  levels),  we  will 
indeed  expect  to  see  how  these  isobars  group  around  certain  points  at  which 
the  isobar  level  becomes  a  local  maximum.  These  points  define  the  locations 
of  the  internal  field's  caustics.  The  higher  the  value  of  a/A  considered, 
the  more  caustics  one  is  expected  to  see,  and  the  more  detailed  and  complex 
the  isobaric  contours  of  the  internal  sound  field  will  be.  These  internal 
isobaric  plots  can  be  constructed  at  increasingly  higher  a/A-  values  to 
observe  the  quantitative  development  of  the  internal  field,  as  well  as  the 
caustic  migrations  for  three  cases:  (1)  when  the  surface  load  is  acting  on  the 
insonified  shell,  (2)  when  it  is  not  acting  on  the  still  insonified  shell,  and 
(3)  when  the  surface  load  is  present,  but  the  shell  is  not  being  insonified  by 
any  incident  wave.  The  interplay  between  these  cases  will  serve  to  derive  a 
variety  of  conclusions,  most  of  them  originally  conceived  on  the  basis  of 
intuition  and  of  our  expectations  for  the  shell  behavior  as  a  mechanical 
filter  of  acoustical  energy. 
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SECTION  4 
NUMERICAL  RESULTS 


The  shell  material  is  initially  steel  of  properties:  (1)  density 

p2  =  7.84  g/cm3,  (2)  dilatational  wave-speed  =  5.8  x  105  cm/sec  and, 

(3)  shear  wave-speed  c„  =  3.1  x  105  cm/sec.  The  shell  is  thin  and  the  a/b 

s2 

ratio  is  initially  180/178  =  1.011235,  where  _a  is  the  outer  shell  radius  and  Jj 
its  inner  radius.  The  shell  is  immersed  in  water  at  20°C,  of  properties:  c^  = 

C  3 

1.476  x  103  cm/sec,  and  pj=  1.0  g/cm  .  The  inner  fluid  is  air  of 
properties:  C3  =  0.344  x  10^  cm/sec,  and  P3  »  0.0012  g/cm^.  The  water  is 
medium  1,  the  shell  is  medium  2,  and  the  air  is  medium  3. 

We  will  consider  two  loads  acting  on  the  outer  shell  surface  r  *  a.  They 
are  shown  in  Figures  1  and  2.  Load  #1,  shown  In  Figure  1,  is  a  wavy  load  that 
acts  on  the  entire  surface  of  the  shell  and  is  described  by  the  relation: 

PL  (8)  *  K  (1  +  -p  cosse) ,  (33) 

where  r  and  s  are  numerical  parameters  associated  with  the  number  of  lobes  in 
the  pattern,  and  the  maxima  or  minima  levels  of  its  amplitude  K.  In  the  above 
text  we  denoted  the  mean  amplitude  K  by  p£.  The  numerical  values  of  r  and  s 
for  this  load  are  r  =  4  and  s  =-  6.  Load  #2,  shown  in  Figure  2,  is  a 
sinusoidal  type  of  load  that  acts  on  a  relatively  small  portion  of  the  shell's 
surface  and  has  equation 
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PL  (0) 


(  0 

<  sin 

l  0 


0  ...  0  <  0  <  ir/3 

(20  -  ■=*■  ) ...  tt/3  <  0  <  5tt/6  . 
...  5ir/6  <  0  <  n 


The  amplitude  of  this  load  was  said  to  be  p2  in  the  text,  but  here  it  is  shown 


to  be  unity  for  reasons  discussed  before.  Both  these  loads  are  symmetric  with 


respect  to  the  0=0  direction,  and  admit  Fourier  cosine  decompositions.  The 


incident  wave  always  hits  these  loaded  shells  from  the  left,  along  the 


0=0  direction. 


We  now  consider  the  Fourier  integral  appearing  in  Equation  (24),  viz. 


T  *  —  f  r  P.  ( ®)  cosn0d0 
n  n  }  o 


and  evaluate  it  for  the  two  loads  given  in  Equations  (33)  and  (34)  and  shown 


in  Figures  1  and  2  .  For  load  #  1  the  result  is. 


P2/P2»  •••  n  =  0 


P2  l 

I_  *  I  (1  +  -r  cos60)  cosnede  =  p2/2rp0,...  n  *  (s=)  6  (36) 

n  "Po  '  0 

0,  ...  any  other  n. 


Thus,  for  r  ■  4,  and  s  =  >6  the  integral  Tn  has  value  P^/Pq  for  n  =  0,  P2/8p0 


for  n  =  6,  and  vanishes  for  any  other  value  of  n.  Clearly,  this  is  the  type 


of  load  that  will  affect  only  a  couple  of  the  shell's  vibrational  modes.  Note 


that  here,  P2/P0  can  be  considered  of  unit  value  since  the  factor  (P2/P0)  was 
already  taken  out  of  the  pertinent  integrals,  as  shown  in  Equation  (25).  For 


load  I  2  we  must  evaluate. 


—  (— )  f5//6  sin  (20  -  2n  /3)  cosn0d0 

If  p  J  IT/  3 


and  the  result  is 
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2  P2 

- 2 - 

n(4-n  )p0 

/T  P2 
&  Po 


Ceos  Ijn  +  COS  ] 


V"  *  2 


for  n  =  2 


This  is  the  type  of  load  that  will  affect  all  the  shell's  modes  of 

A 

vibration.  The  values  of  In,  as  given  in  Equation  (27),  are  merely  the  above 

2  2 

values  of  7  in  Equations  (36)  and  (38),  times  the  factor  -  (P1/P2)  kc  a  . 

All  the  coefficients  have  now  been  evaluated,  and  we  can  proceed  to  display 
computed  plots.  We  have  developed  computer  codes®’®  that  permit  the 
evaluation  of  unloaded  shell  responses,®"*®  viewed  in  the  outer  medium  1.  We 
extended  these  codes  to  also  compute  the  internally  transmitted  pressure 
fields  into  loaded  shells  as  we  have  described  it  here.  Additional  codes  were 
created  to  compute  and  graph  the  isobars  (i.e.,  the  loci  of  points  exhibiting 
the  same  value  of  internal  pressure)  inside  the  shell,  in  various  cases,  and 
as  function  of  all  the  problem  parameters. 

Figure  3  shows  four  successive  stages  of  development  of  the  internal 
isobaric  fields  inside  the  shell  for  arbitrary  values  of  k3a  of  1,  3,  5,  and 
10.  These  patterns  are  computed  for  an  unloaded  steel  shell  using  Equation 
(14')  and  (14).  A  point  (r,  6)  defines  any  location  in  medium  3,  inside  the 
shell.  The  allowable  ranges  are  0  <  r  <  b  and  0  <  8  <  n  .  The  patterns  are 
symmetric  with  respect  to  the  centerline  direction,  e  =  0.  The  incident  wave 
of  amplitude  pg  always  hits  the  shell  from  the  left.  As  we  anticipated,  the 
internal  fields  become  more  and  more  developed  and  complex  as  the  k3a-value  is 
increased.  (Note:  k3  =  u>/C3  =  2irf/C3,  and  c3  is  the  sound  speed  in  air). 

For  k3a  =  10  the  pattern  has  already  about  eight  noticeable  caustics,  while  at 
k3a  *  5,  there  are  two  noticeable  ones  on  the  centerline,  at  polar 
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coordinates  of  (0.3  ,  0°)  and  (0.4,  180°).  The  radial  coordinates  have  been 
nondimensionalized  to  r/b.  The  artificial  irregularities  in  the  computed 
patterns  near  their  outer  edges  are  numerically  caused  by  the  size  of  the  step 
that  separates  consecutive  isobars  within  each  plot. 

Figure  4  shows  the  same  four  isobaric  contour  plots  inside  the  same 
shell,  at  the  same  values  k3a  =  1 ,  3,  5  and  10,  now  computed  for  a  shell 
loaded  with  Load  #1  (see  Figure  1).  The  lower-frequency  plots  (i.e.,  k3a  =  1 
and  3)  are  not  too  informative  since  hardly  any  caustic  has  had  a  chance  to 
develop  yet,  but  for  k3a  =  5,  and  10  we  already  see  a  fully  developed  isobaric 
set  of  contours  not  very  different  from  those  in  Figure  3  for  an  unloaded 
shell  at  the  same  k3a  -values.  For  k3a  =  5,  the  same  two  caustics  are 
observed  at  the  centerline,  at  about  the  same  places.  The  similarity  between 
these  plots  and  those  for  the  unloaded  shell  in  Figure  3  is  undoubtedly  due  to 
the  fact  that  Load  #1  is  a  type  of  load  that  affects  only  a  couple  of  the 
shell's  vibrational  modes  and  thus,  it  produces  no  overwhelming  distortion  of 
the  pattern  already  present  in  the  unloaded  shell  case.  The  plots  are 
computed  from  Equation  (21)  with  qn  as  in  Equation  (18),  and  the  integral 
within  coefficient  %i*  as  in  Equation  (36);  in  other  words. 


...for  n  =  0 

...for  n=6  (39) 

...for  any  other  n. 


and  where,  in  addition,  we  have  assumed  p£/pg  =  1  to  generate  the  plots.  This 
assumption  was  introduced  to  show  that  even  when  the  amplitude  of  the  surface 
load  P2  is  comparable  to  that  of  the  incident  wave  pg,  its  effect  on  the 
internal  pressure  field  generated  by  the  amplitude  pQ  of  the  incident  wave,  is 
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minor.  This  is  due  to  the  peculiar  symmetric  nature  of  the  surface  load 
distribution  of  Load  #  1. 

Figure  5  exhibits  the  isobaric  contours  inside  the  shell,  at  non- 
dimensional  frequencies  of  k3a  *  3,  5,  and  10,  but  now  for  a  shell  loaded  with 
Load  #  2  (see  Figure  2).  Fully  developed  internal  isobaric  fields  start  to  be 
seen  for  k3a  =  5  and  more  so  for  k3a  =  10.  Even  though  Load  #  2  is  a  type  of 
load  that  affects  all  the  shell's  vibrational  modes,  the  caustic  distribution 
inside  the  shell  is  quite  similar  to  those  of  the  unloaded  shell  (c/o  Figure 
3)  at  the  same  frequencies.  The  distortions  are  now  greater  than  for  Load  #  1 
(Figure  4),  but  the  caustics  remain  more  or  less  in  the  same  places,  although 
their  local  value  changes  somewhat.  This  is  quite  obvious  at  k3a  =  10,  and 
less  evident  for  k3a  =  5  where  only  the  right  caustic  on  the  centerline 
remains  in  place,  while  the  other  seems  to  have  migrated  toward  the  left -most 
edge  of  the  shell.  This  figure  is  computed  just  as  Figure  4,  but  the 
coefficient  %\*  now  is  given  by  the  expression 


At 


—  kc  a  { J -j ( k i a )  +  (■—-)  - z  Ceos  +  cos  ]} 

P2  s2  1  nv  1  W  b  3  J 

Pi  2  2  P2  Jf 

—  kSj  a  CJn(k,a)  ♦  (pj>  <-  •*!  >  >  *  f°r  "  "  2 


Vn*2 

(40) 


We  also  chose  P2/P0  *  1  To  generate  Figure  6  as  we  did  in  Figure  4. 

It  seems  that  k3a  =  5  is  the  (nondimensional )  frequency  value  considered 
here  at  which  the  cleanest  and  easiest-to-interpret  results  are  found.  Higher 
k3a  -values  clearly  clutter  the  internal  pattern  too  much  for  the  extraction 
of  simple  conclusions.  The  generation  of  many  of  these  graphs  can  be  quite 
expensive. 
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We  still  have  to  consider  the  effect  of  the  surface  load  al one  on  the 
internally  transmitted  pressure  field.  This  can  be  done  by  the  superposition 
principle  and  the  results  derived  in  Section  3.  We  still  use  P2/P0  =  1  to 
generate  the  results,  which  are  obtained  by  means  of  Equation  (30),  with  qn  as 
given  in  Equation  (26).  The  integral  I  is  given  by 


1 

n 


Pl  2 
—  k 
P  2  s2 


2 


a 


(41) 


where  I  is  given  in  Equation  (36)  for  Load  #1 ,  and  in  Equation  (38)  for  Load 
#2.  About  20  terms  were  added  in  the  series  of  Equation  (30).  Figure  6  shows 
the  result  of  this  calculation  for  the  single  value  of  k3a  =  5  for  both  Load 
#1  (upper  plot)  and  for  Load  #2  (bottom  plot).  Since  only  the  load  is  acting 
(there  is  no  'ncident  wave  now),  the  upper  plot  is  symmetric  about  the  origin, 
since  it  is  caused  by  the  axially  symmetric  Load  #1.  Four  caustics  are 
observed,  one  fully  developed  at  the  origin,  and  three  more  starting  to  grow 
around  the  shell's  inner  surface.  The  .ower  plot  is  asymmetric,  since  it  is 
generated  by  the  non-axially  symmetric  Load  #2.  The  caustic  at  the  origin, 
observable  in  the  upper  plot,  has  spread  out  to  a  wider  region  in  the  bottom 
plot. 

We  can  now  quantitatively  address  the  question  of  the  advantage  of 
sensing  transmitted  fields  inside  the  shell,  in  contrast  to  fields  sensed 
outside  it.  Equation  (31)  gave  the  theoretical  ratio  PR  for  this  advantage. 

The  following  discussion  will  be  limited  to  the  case  kja  =  5,  for  which 
all  the  ingredients  have  already  been  determined  for  both  types  of  loads.  The 
unloaded  shell  case  and  its  isobaric  contours  were  displayed  in  Figure  3.  The 
denominator  corresponds  to  the  case  of  a  surface-loaded  shell  without  an 
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incident  wave,  a  situation  displayed  in  Figure  6  for  Load  #1  (upper  plot)  and 
for  Load  #2  (bottom  graph).  It  does  not  make  much  sense  to  compute  values  for 
the  ratio  at  every  point  inside  the  shell.  It  is  more  pertinent  to  consider 
the  values  of  the  PR  -ratio  only  at  the  locations  of  the  caustics  of  its 
numerator  (viz,  those  shown  in  Figure  3).  The  two  caustics  of  Figure  3  are 
located  at  (0.3,  0°)  and  (0.4,  180°)  and  have  values  of  about  1.3  and  1.2, 
respectively.  For  Load  #1  the  upper  graph  in  Figure  6  shows  values  of  about 
0.18  and  0.38,  respectively,  at  the  above  mentioned  caustic  locations.  For 
load  #2  the  lower  graph  of  Figure  6  shows  values  of  0.50  and  0.74, 
respectively  at  those  two  caustics.  The  advantage  ratio  PR  in  all  cases  is 
summarized  in  Table  1. 

TABLE  1.  VALUES  OF  PR  [EQUATION  (31)1  for  k3a  -  5 


Caustic  # 

Load  # 

Left  Caustic  (0.4,  180°) 

Right  Caustic  (0.3,  0°) 

Load  #  1 

=  6.66  (16.5  dB) 

=  3.42  (10.7  dB) 

Load  #2 

jffa  =1.62  (4.2  dB) 

^|Q=  2.60  (8.3  dB) 

The  values  in  decibels,  shown  in  parenthesis  in  Table  1,  are  20  log  of  the 
value  of  the  shown  ratio  in  each  case.  It  is  clear  that  the  ratios  are  larger 
than  one  in  all  cases,  and  depending  on  which  load  and  which  caustic  one 
considers.  Table  1  shows  advantages  ranging  from  4  to  about  16  dB.  For  the 
construction  of  the  values  in  Table  1,  we  assumed  that  pg/p?  =  1.  Had  this 
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external  "signal -to-noise"  ratio  been  smaller,  the  advantages  in  Table  1  would 
have  been  even  larger  than  shown  there. 

We  can  analyze  the  internal  SPL  as  a  function  of  frequency,  at  selected 
points  inside  the  shell,  for  the  various  loaded  or  unloaded  conditions  we 
considered.  We  could  also  investigate  the  SPL  as  a 

function  of  radial  location  across  a  shell  diameter,  at  selected  values  of  the 
frequency,  for  all  the  loading/unloading  conditions. 

Figures  7a  and  7b  show  the  variation  of  the  (normalized)  internal  SPL: 
(i.e.,  P3/P0)  as  a  function  of  k3a  in  the  range  0  <  k3a  <  10,  at  11  points  in 
the  upper  half  of  the  shell's  interior.  If  (r,0)  are  plane  polar  coordinates, 
these  11  points  are,  for  Figures  (7a):  (b,0°),  (b,45°),  (b,90°),  (b,135°)  and 
(b,  180°)  and  for  Figure  (7b):  (0,0°),  (b/2,  0°),  (b/2,  45°),  (b/2,  90°), 

(b/2,  135°)  and  (b/2,  180°).  These  results  are  for  a  cylindrical  unloaded 
shell  and  they  are  evaluated  by  means  of  Equations  (14')  and  (14).  The 
significance  of  all  these  peaks  can  be  given  in  terms  of  shell  resonances. 
These  resonances  could  be  investigated,  one  mode  at  the  time,  at  a  later 
date.  As  k3a  is  increased,  the  normalized  SPL,  | P3/P0 I »  1S  basically  zero 
until  a  shell  resonance  is  excited  and  a  peak  is  generated.  These  peaks  are 
seen  to  be  very  steep  and  at  any  one  of  them  the  value  of  P3  can  be  ten  or 
even  one  hundred  times  larger  than  the  incident  wave  amplitude  pg.  Thus,  it 
follows  that  sound  transmission  into  the  shell  can  only  occur  through  narrow 
spectral  windows  centered  at  the  shell's  modal  resonances,  and  that  the 
amplitude  of  the  transmitted  fields  can  be  so  large  that  the  shell  can  behave 
as  a  very  effective  sound  concentrator.  It  is  not  simple  to  analytically 
extract  these  conclusions  from  all  the  (11)  points  shown  in  Figure  7,  but  it 
is  relatively  simple  to  show  them  for  the  point  at  the  center  of  the  shell 
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(r=0),  where  the  only  two  peaks  that  are  present  occur  at  values  of  kga  =  3,8 
and  7.0,  which  are  the  first  two  roots  of  Jg'^a)  =  0  in  the  considered 
range.  This  will  be  demonstrated  toward  the  end  of  this  section. 

Figures  8a  and  8b  show  the  same  frequency-dependent  plots  evaluated  at 
the  same  11  points  inside  the  shell  but  now  for  a  shell  loaded  with  Load 
#1 .  The  shell  is  excited  by  both  surface  load  and  incident  wave.  The 
coordinates  of  the  11  interior  points  are  given  above  and  the  pertinent 
formulas  used  in  the  calculation  are  given  in  Equations  (32),  (14),  and 
(26).  We  assumed  that  P2/P0  =  1  To  generate  these  plots.  The  coefficients 
I  appearing  in  Equation  (26)  are  given  in  Equations  (41)  and  (36). 
Comparisons  of  the  corresponding  portions  of  Figures  7  and  8  show  that  in  both 
cases  the  resonance  peaks  appear  approximately  at  the  same  locations  with 
small  changes  in  amplitude  and  with  the  occasional  emergence  of  an  additional 
peak.  All  these  differences  are  minor  and  can  be  attributed  to  the  central 
symmetry  of  Load  #1,  which  we  saw  was  of  the  type  affecting  only  a  couple  of 
the  shell's  modes.  Therefore,  it  is  not  expected  that  these  will  change  very 
drastically  the  response  of  the  unloaded  shell  to  the  incident  wave.  We  now 
understand  better  the  reason  for  the  similarity  between  Figures  7  and  8. 

Figures  9a  and  9b  exhibit  the  same  plots  shown  in  Figure  8,  but  now  they 
are  computed  for  the  case  where  the  shell  is  loaded  with  Load  #2.  The 
formulas  required  to  generate  these  plots  are  the  same  as  for  Figure  8,  but 
now  the  load-dependent  coefficients  I  are  given  by  Equations  (41)  and 
(38).  Figure  9  is  computed  at  the  same  eleven  internal  points  as  Figures  7 
and  8.  The  two  resonance  spikes  at  the  center  of  the  shell  (r  =  0)  mentioned 
in  Figure  7,  remain  at  the  same  places  but  their  amplitudes  seem  to  have 
changed.  For  most  of  the  other  ten  points,  the  spectra  shown  in  Figure  9 
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contain  more  (and  different)  resonance  peaks  than  those  shown  in  Figure  7  for 
the  unloaded  shell  case.  Load  #2  is  the  type  of  load  that  affects  all  the 
shell  modes,  and  it  is  not  surprising  that  the  same  ten  internal  points  now 
show  more  complicated  pressure  spectra,  containing  more  resonance  peaks  than 
before.  These  additional  peaks  are  clearly  due  to  the  nature  of  the  surface 
load.  The  information  in  these  plots  already  permits  the  identification  of 
which  of  the  resonance  peaks  at  which  of  the  internal  points  is  caused  by  the 
surface  load,  and  which  ones  by  the  incident  wave.  Future,  more  detailed 
studies  will  allow  us  to  evaluate  the  distorting  effect  of  the  surface  load  on 
the  shell's  modal  resonances  as  excited  by  the  incident  wave,  for  values  of 
the  external  "signal -to-noise"  ratio,  Pq/P2»  different  than  unity.  This  ratio 
was  also  assumed  unity  in  the  figures  that  follow. 

Figure  10  begins  a  study  of  the  monopole  (i.e.,  n  =  0)  mode  of  vibration 
of  the  shell  when  it  is  unloaded.  For  this  symmetrically  breathing  mode,  we 
plot  the  contribution  of  the  n  =  0  partial -wave  or  normal  mode  to  the 
(normalized)  internal  SPL,  | P3/P0 | »  versus  kga,  at  six  radial  locations  (r  = 

0,  0.2b,  0.4b,  0.6b,  0,8b,  and  b)  along  any  radius.  The  result  is  the  same  in 
all  directions  because  for  n  =  0,  cosne  =  0,  for  all  6.  This  is  a  plot  of 
the  first  (i.e.,  the  n  =  0)  term  in  the  series  in  Equation  (14').  The  two 
peaks  already  observed  before  at  kja  =3.8  and  7.0  in  Figure  7  are  present  at 
all  radial  locations  with  an  amplitude  somewhat  different  for  each  radial 
location.  At  r  =  0,  the  two  peak  amplitudes  are  the  same  as  those  already 
seen  in  Figure  7b,  and  for  other  locations  they  change  as  shown.  All  the 
peaks  at  kja  =3.8  and  7.0  seen  in  all  the  plots  of  Figure  7a  (for  r  =  b)  are 
caused  by  this  monopole  (n  =  0)  mode,  as  can  be  easily  seen  in  the  top  graph 
of  Figure  10. 
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Figure  11  is  the  counterpart  of  Figure  10,  but  now  for  the  case  of  a 
shell  loaded  with  Load  #1.  Figure  11  displays  plots  of  the  first 
(i.e.,  n  =  0)  term  of  Equation  (32)  versus  k3a.  Here  qn,  qn  are  given  by 

A 

Equations  (14)  and  (26)  and  where  the  only  remaining  In  coefficient  is 
-  2  2 

I0  *  (“P1/P2)  k  a  •  The  same  two  resonance  peaks  present  in  Figure  10  still 
s2 

persist  in  Figure  11  but  with  quite  different  amplitudes  than  before.  We 

recall  that  Load  #1  affects  only  a  couple  of  the  modes,  but  the  monopole  mode 

(n  =  0)  is  one  of  them,  and  this  explains  the  substantial  amplitude  difference 

seen  between  the  resonance  peaks  of  Figures  10  and  11. 

Figure  12  displays  the  same  results  as  Figure  11  but  now  for  the  shell 

loaded  with  the  substantially  different  Load  #2.  The  formulas  used  in  the 

calculations  are  as  in  the  above  paragraph  except  that  now 
2  2 

lo  =  (-P1/P2)  a  /*•  We  note  that  the  same  persistent  peaks  appearing 
s2 

before  at  k3a  =  3.8  and  7.0  continue  to  be  present  for  all  radii  but  now  with 
amplitudes  quite  different  from  those  in  Figure  10,  and  of  course,  from  those 
in  Figure  11. 

The  effect  of  the  monopole  resonances  is  very  persistent  and  strong.^ 

All  the  computed  plots  of  the  summed  interior  SPL  consistently  exhibit  those 
two  strong  narrow  peaks  at  the  mentioned  pair  of  spectral  locations.  Peaks  at 
other  frequencies  are  caused  by  other  higher  modes.  At  values  of  k3a  away 
from  those  resonance  values,  the  (otherwise  strong)  monopole  contribution  to 
the  fields,  disappears  and  the  fields  become  mixtures  of  (weaker)  higher-order 
modes.  We  see  now,  for  example,  that  all  the  isobaric  contour  plots  shown  in 
Figures  3  to  6  and  computed  for  arbitrarily  chosen  values  of  the  non 
dimensional  frequency  of  k3a  =  1,  3,  5,  and  10,  have  all  missed  the  monopole 
resonance  contributions  which  occur  within  narrow  bands  around  3.8  and  7.0. 
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This  merely  means  that  the  internal  fields  and  values  of  the  advantage  ratio 
Pr  computed  earlier,  could  be  considerably  enhanced  if  they  were  recomputed  at 
the  preferred  k3a-values  of  3,8  and  7.0.  Therefore,  the  natural  filter  and 
lens-like  or  sound  concentrating  effect  of  the  shell  can  be  enhanced  beyond 
the  values  shown  earlier  (i.e.,  in  Table  1),  by  operating  near  these  two 
resonances  of  the  monopole  mode  contained  within  the  original  broad  band 
(0  <  k3a  <  10)  considered. 

The  monopole  mode  contribution  to  the  internal  fields  can  also  be  studied 

by  computing  its  radial  variation  for  various  values  of  k3a.  Figures  13  to  15 

summarize  this  effect  and  its  radial  dependence  just  as  Figures  10  to  12  studied 

it  for  its  frequency  dependence.  Figure  13  is  a  plot  of  the  internal  pressure 

associated  with  the  first  (i.e.,  n  =  0)  term  of  the  series  in  Equation  (14*) 

for  an  unloaded  shell  displayed  versus  radial  coordinate  r  in  the  range 

0  <  r/b  <  1 ,  at  the  usual  four  arbitrary  values  of  k3a  (viz.,  k3a  =  1,  3,  5 

and  10)  that  we  have  chosen  to  span  the  band  0  <  k3a  <  10.  Figure  14  is 

analogous  to  13  and  it  exhibits  the  first  (i.e.,  n  =  0)  term  of  Equation  (32), 

-  ,  2  2 

for  Load  #1.  The  coefficient  In  here  is  (-P1/P2)  k_  a  as  we  found  for 

s2 

Figure  11,  and  we  have  assumed  P2/P0  =1.  Figure  15  displays  the  first  term 

of  the  series  in  Equation  (32),  now  for  Load  #2.  In  this  case  the  coefficient 
-  2  2 

Io  is  (-P1/P2)  k  a  /n  as  it  was  for  Figure  12,  and  we  still  assume  p2/Po  = 

1.  When  we  compare  the  plots  in  these  three  figures  for  corresponding  values 
of  k3a,  we  see  that  the  transmitted  sound  pressure  amplitudes  vary  quite 
substantially  as  we  go  from  Figure  13  (unloaded  case),  to  Figure  14  (Load  #1) 
to  finally  Figure  15  (Load  #2).  Both  these  loads  affect  the  monopole  (n  =  0) 
mode  and  this  explains  the  variations  in  corresponding  internal  pressure 
amplitudes  observed  between  these  three  figures.  It  is  quite  noticeable  that 
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the  nulls  present  in  these  plots  at  each  k3a-value  do  not  change  their  radial 
locations  as  we  go  from  Figures  13  to  15.  Closer  inspection  shows  that  in  the 
range  of  the  figures  (i.e.,  0  <  k3a  <  10),  the  nulls  seen  in  Figure  13  are  all 
manifestations  of  the  first  three  zeroes  of  Jg(k3r),  which  are  known  to  occur 
at  k3r  =  2.40,  5.52,  and  8.65.  Since  k3a  =  10  for,  say,  the  top  graph  in 
Figure  13,  division  of  these  zeroes  by  ten  yields  r/a-values  of  0.24,  0.55  and 
0.86.  These  r/a-values  are  undistinguishable  from  the  r/b-values  shown  as  the 
abscissas  of  the  nulls  in  that  plot  because  the  a/b-value  for  this  shell  is 
very  close  to  unity  (viz,  1.0112).  For  the  second  plot  in  Figure  13, 
calculated  for  k3a  =  5,  division  of  the  three  roots  by  five  eventually  yields: 
r/b  =  0.48,  1.10,  and  1.73.  Only  the  first  of  these  appears  in  the  range 
(i.e.,  0  <  r/b  <  1)  where  the  plot  was  computed.  Analogous  results  are  found 
for  all  other  plots  of  Figure  13  and  for  those  in  Figures  14  and  15.  It 
follows  that  these  nulls  in  the  monopole  contribution  to  the  internal  field 
are  controlled  by  the  zeroes  of  Jg  (k3r)  since  this  factor  is  indeed  all  that 
remains  in  Equation  (14* ),  for  the  monopole  (i.e.,  n  =  0)  case. 

It  remains  to  formally  demonstrate  that  at  the  center  of  the  shell  (i.e., 
r  =  0)  the  resonance  peaks  appearing  in  the  spectral  plots  are  controlled  by 
the  resonance  conditions  Jg'(k3a)  =  0.  Since  the  6x6  determinants  involved 
here  were  so  complicated,  it  was  unexpected  to  find  such  a  simple  condition. 
This  result  was  already  anticipated  when  we  discussed  Figure  7. 

For  r  =0  all  the  terms  in  the  series  shown  in  Equation  (14')  vanish 
except  for  the  n  =  0  term.  This  is  so  because  all  the  Bessel  functions  of 
zero  argument  vanish  except  Jg(0)  which  takes  on  value  unity.  The  n  *  0  term 
reduces  to  coefficient  qg.  Using  the  entries  in  Appendix  A  this  coefficient 
takes  the  form 


(42) 


Pi  2  2  0  0 

—  k  a  Jo  (k1a)Oig  +  ^l3  Jo,(*cia)02B 

q  =  - 

O  P3  2  2  0  0 

—  k  b  J0  (k3b)0,6  +  k  3b  Jo'(k3b)D56 

o  o 

where  Di6,  D26.  ...  are  the  5x5  minor  determinants  of  the  elements  d16,  d^g, 
...  listed  in  Appendix  A,  all  evaluated  for  n  =  0.  Resonances  in  the  SPL  are 
associated  with  the  vanishing  of  the  denominators  of  the  series.  For  this 
single  term  in  Equation  (42),  resonances  are  associated  with  the  vanishing  of 
its  denominator.  Since  the  inner  fluid  (air)  has  a  density  p3  very  tenuous 
compared  to  that  of  the  steel  shell  (i.e.,  p3  <<  P2),  it  follows  that  unless 
the  shell  is  extremely  thin,  the  first  term  in  the  denominator  of  Equation 
(42)  is  negligible  compared  to  the  second,  and  that  second  one  immediately 
yields  the  condition  Jo'(k3b)  =  0.  Since  the  b/a-ratio  for  this  shell  is  very 
close  to  unity,  this  condition  is  undistinguishable  from  JQ'  (k3a)  =  0,  which 
yields  the  results  observed  in  Figure  7,  and  which  we  wanted  to  formally 
demonstrate  here  for  this  air-filled  thin  steel  shell. 

Basic  explanations  for  the  origin  of  the  internal  caustics  at  the  higher 
k3a  -  values,  for  which  the  internally  transmitted  fields  are  more  developed, 
can  be  given  in  terms  of  creeping  waves  as  we  did  earlier  for  coated 
cylindrical  shells,6  and  for  fluid-filled  cavities  within  elastic  media. H 
is  well  known  that  transmission  of  sound  through  (curved)  shells  differs 
significantly  from  sound  transmission  through  flat  plates  only  at  the  low-end 
of  the  frequency  spectrum  (i.e.,  in  the  Rayleigh  region  and  in  a  portion  of 
the  resonance  region).  The  previously  mentioned  cut-off  value  of  about  k3a- 
20,  where  the  numerical  instabilities  and  the  slowness  of  the  series 
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convergence  begin  to  prevent  the  possibilities  of  computation,  seems  to  be  the 
natural  upper  bound  for  the  spectral  region  beyond  which  sound  transmission 
through  shells  could  well  be  modeled  by  means  of  flat  plate  theories.  All  our 
sound  structure  interaction  work  that  we  have  quoted  above  (as  well  as  other 
references  not  mentioned  here)  was  summarized  in  another  (updated)  review  of 
the  Resonance  Scattering  Theory  (RST)  that  appeared  in  1982.^ 

We  conclude  this  sub-section  with  a  statement  of  the  version  of  the 
Fourier  theorem  that  permits  the  decomposition  of  the  arbitrary  load  acting  on 
the  shell  given  in  Equation  (16). 

Theorem:  Let  f(e)  be  a  function  satisfying  Dirichlet  boundary  conditions  on 
0  <  0  <  a,  and  let 

?c  (n)  =  Jo  f(B)  cos-~de  (43a) 

be  its  Fourier  finite-cosine  transform.  Then  f(0)  admits  the  representation 

f(0)  *7  I  en  f.  (n)  cos (43b) 

a  n=0  n  c  a 

at  each  point  of  (0,a)  where  f  (0)  is  continuous. 

Corolary :  Let  a  =  if  ,  then  at  each  point  of  (0,  v)  the  above  finite 
transform  pair  is 


f  (n)  =  / o  f( 0)  cosnode 


(44a) 


f(e) 


i 

IT 


Jo  en  fc(n)  cosne- 


(44b) 


31/32 


NSWC  TR  85-42 


► 

i 


j* 


K 


SECTION  5 

SHELL  THICKNESS  AND  MATERIAL  COMPOSITION  EFFECTS 


This  section  presents  a  detailed  study  of  the  effects  of  shell  thickness 
and  material  composition  on  the  pressure  fields  transmitted  inside  the  doubly 
excited  shell.  The  computed  displays  presented  up  to  this  point  pertain  only 
to  thin  shells  of  a  single  thickness  and  composition.  Those  results  will  now 
be  extended  to  include  thicknesses  as  high  as  20  percent,  and  compositions 
including  aluminum  and  steel.  Roth  an  incident  plane  wave  and  an  external 
surface  excitation  act  on  the  outer  surface  of  the  submerged  shell. 

For  the  reasons  described  earlier,  we  will  select  a  fixed  nondimensional 
frequency  such  that  k3a*5.0.  For  this  frequency,  it  follows  that  X«a,  which 
means  that  one  is  at  the  edge  of  the  resonance  region.  For  any  other  selected 
frequency  the  calculations  follow  a  similar  pattern,  although  at  high 
frequencies,  both  analysis  and  computations  are  harder  to  carry  on.  It  is  of 
interest  to  first  study  the  isobaric  contour  plots  of  the  insonified  shell  in 
the  absence  of  surface  loads  and  later  to  recompute  them  when  the  submerged 
shell  is  additionally  excited  by  Loads  #1  and  #2. 

Figure  16  shows  the  isobaric  contours  of  the  internal  SPL  when  the 
submerged  shell  is  insonified  but  subject  to  no  surface  excitation.  The  shell 
material  is  steel,  the  frequency  is  such  that  k3a=5.0,  and  the  thicknesses  are 
such  that  b/a=0.80,  0.85,  0.90  and  0.95.  The  calculations  are  performed  by 
means  of  Equations  (14')  and  (14).  The  thicker  the  shell  is,  the  smaller  the 
local  maxima  that  are  observed  at  the  various  corresponding  caustic 

33 


locations.  This  merely  means  that  away  from  resonances,  thicker  shells  allow 
for  the  transmission  of  smaller  amounts  of  sound  energy  into  its  internal 
points.  The  caustic  locations  shift  with  thickness  but  in  a  slow  fashion  that 
easily  permits  the  tracing  of  a  given  caustic  as  the  shell  thickness  is 
changed. 

Figure  17  shows  analogous  contour  plots  for  a  steel -shell  subjected  to 
the  incident  pressure  wave  and  to  the  surface  excitation  generated  by  Load 
#1.  The  nondimensional  frequency  is  still  k3a=5.0,  and  the  same  four  shell 
thicknesses  (i.e.,  b/a=0.80,  0.85,  0.90,  and  0.95)  are  displayed.  These  plots 
are  not  too  different  from  the  corresponding  ones  in  Figure  16,  particularly 
for  the  thinner  shells.  These  contours  are  computed  using  Equations  (21), 
(18),  (20),  and  (36).  The  similarity  between  corresponding  plots  in  Figures 
16  and  17  can  be  explained  because,  as  we  have  seen.  Load  #1  merely  affects  a 
couple  of  the  shell  modes. 

Figure  18  displays  the  same  information  contained  in  Figure  17  but  now 
the  outer  shell's  surface  is  excited  with  Load  #2.  The  differences  between 
Figures  18  and  16  are  more  pronounced  than  between  Figures  17  and  16  because 
Load  #2  affects  all  the  shell  modes. 

The  information  exhibited  in  Figures  16  to  18  for  k3a=5  is  given  in 
Figures  19  to  21  for  k3a=7.0.  In  Figure  19  the  shell  is  insonified  but  there 
is  no  surface  load;  in  Figure  20  the  shell  is  excited  with  Load  #1;  and  in 
Figure  21  it  is  excited  with  Load  #2.  At  higher  frequencies  one  observes  more 
caustics  that  are  more  developed  than  at  lower  frequencies. 

Steel  and  aluminum  do  not  behave  in  the  same  fashion,  and  the  differences 
are  displayed  in  Figures  22  to  27.  These  six  figures  display  information 
about  the  internally  transmitted  SPL  analogous  to  that  shown  in  Figures  16 


to  21,  but  now  for  an  aluminum  shell.  Figure  22  shows  the  internal  SPL 
contours  for  the  usual  four  thicknesses,  at  the  fixed  nondimensional  frequency 
of  kja^B  for  an  unloaded  aluminum  shell.  The  aluminum  properties  used  in 
these  plots  are  :  =2.70  g/cm,3  0^^=6.374x10^  cm/sec,  and  c^S.lllxlO5 

cm/ sec.  The  calculations  are  performed  using  Equations  (14)  and  ( 14 ' ) . 

Figure  23  displays  similar  results  when  the  insonified  shell  is  excited  by  the 
outer  surface  Load  #1.  In  this  case,  the  calculations  are  done  by  means  of 
Equations  (21),  (18),  (20),  and  (36).  Figure  24  exhibits  the  corresponding 
results  for  the  case  of  Load  #2,  in  which  case,  the  calculations  are  made  by 
means  of  Equations  (21),  (18),  (20),  and  (38).  Figure  25  deals  with  the 
insonified  but  otherwise  unloaded  aluminum  shell  in  water  at  k3a=7.0.  The 
differences  in  transmissibility  due  to  steel  and  aluminum  emerge  from  the 
comparison  of  Figures  19  and  25.  Aluminum  being  a  "softer,"  more  penetrable 
material  allows  the  transmission  of  stronger  internal  SPL's  than  steel,  at 
about  the  same  internal  points,  and  for  corresponding  shell  thicknesses. 

Figure  26  corresponds  to  the  aluminum  shell  at  k3a=7.0  when  the  shell's 
surface  is  excited  with  the  symmetric  Load  #1.  Figure  27  displays  analogous 
information  when  the  Load  #2  is  the  one  exciting  the  outer  surface  of  the 
insonified  aluminum  shell.  Comparing  Figures  27  and  21  for  aluminum  and 
steel,  respectively,  shows  that  for  the  same  thicknesses,  the  stiffer  steel 
shell  permits  the  transmissibil ity  of  a  weaker  internal  pressure  field  than 
aluminum,  at  about  the  same  spatial  points.  Very  lengthy  quantitative 
comparisons  could  be  extracted  from  the  detailed  comparisons  of  all  the  cases 
displayed  in  all  these  figures.  We  have  just  stated  the  most  salient  global 
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features  which  are  most  evidently  noticeable  from  a  quick  inspection  of  the 
displayed  results. 

The  effect  of  the  isolated  surface  excitation  (either  Load  #1  or  #2)  on  a 
steel  or  aluminum  shell,  is  shown  in  Figures  28  to  31.  By  means  of  the 
superposition  principle  and  the  analysis  of  the  doubly  excited  shell  in 
Section  3,  it  is  possible  to  study  the  effect  of  the  surface  load  in  the 
absence  of  the  incident  wave.  The  calculations  are  performed  by  means  of 
Equations  (30),  (26),  and  either  Equation  (36)  for  Load  #1  or  Equation  (38) 
for  Load  #2.  Figure  28  displays  the  internally  transmitted  pressure  contours 
due  to  the  isolated  surface  excitation  of  Load  #1  acting  on  the  steel  shell. 
Figure  29  exhibits  the  same  information  for  solely  Load  #2  acting  on  the  steel 
shell.  Figure  30  pertains  to  Load  #1  acting  alone  on  the  outer  surface  of  the 
aluminum  shell.  Figure  31  exhibits  the  same  information  when  the  shell 
excitation  is  caused  by  Load  #2.  All  these  plots  are  computed  for  the  usual 
four  shell  thicknesses  (viz,  b/a=0.80,  0.85,  0.90,  and  0.95)  at  the  fixed  non- 
dimensional  frequency  of  k3a=5.0.  These  results  are  the  ones  needed  to 
compute  the  sum  appearing  in  the  denominator  of  Equation  (31),  which  yields 
the  advantage  ratio  PR.  A  summary  of  results  for  all  the  thicknesses,  loads, 
and  shell  compositions  is  given  in  Table  2.  This  table  also  shows  the  results 
of  computing  the  advantage  ratio  PR  in  all  cases,  just  as  we  did  it  earlier  in 
Table  1.  The  ratio  PR  almost  always  comes  out  to  be  greater  than  one,  and 
thus,  it  follows  that  for  these  material  compositions  and  shell  thicknesses, 
under  the  shown  conditions  of  insonification  and  external  excitation,  it  still 
appears  to  be  more  advantageous  to  sense  pressure  fields  inside  rather  than 
outside  the  shell.  Ideal  spatial  locations  have  been  identified  by  means  of 
the  caustics  generated  by  the  overall  lense-like  action  of  the  shell. 
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TABLE  2.  VALUES  OF  PR  FOR  VARIOUS  SHELL-THICKNESSES  AND  COMPOSITIONS 

The  advantage  ratio  Pp  for  various  shell  thicknesses  and  compositions.  From 
Eq.  (31),  using  k3a=5.0.  Steel  above,  aluminum,  below. 


STEEL 

b/a= 

0.80 

0.85 

0.90 

0.95 

LOAD  # 


LOAD  1 
LOAD  2 


LOAD  1 
LOAD  2 


LOAD  1 
LOAD  2 


LEFT  CAUSTIC  ON  AXIS 


0.18  _  1.0  (OdB) 

!m 

0.18  1.2  (1.6dB) 

OTIS'  * 


0.275  _  4.6  (13.2dB) 
OTW  " 

0.275  _  2.75  (8.8dB) 
OTTO 


RIGHT 

CAUSTIC  ON  AXIS 

0.22 

1.8  (5.3  dB) 

tr.iz  ~ 

0.22 

1.5  (3.4dB) 

OTIS'  " 

0.30 

7.5  ( 17 . 6dB) 

0T0I~ 

0.30  _ 

0.008 

0.225 

_  4.5  (13. dB) 

OTOSO" 

0.225 

_  2.25  (7. OdB) 

or 

0.50  _ 

7.1  ( 17. ldB) 

0T07 

0.50  _ 

.  2.0  (6. OdB) 

OTZF 

RIGHT 

CAUSTIC 

1  5.0 

40  (32dB) 

urns- 

5.0 

0.2 

0.55 

_  3.7  ( 11. 3dB) 

OTIS' 

0.55 

,  0.07 

ST0~ 

0.75 

„  5.0  ( 14dB) 

om 

0.75 

,1.5  (3.5dB) 

OTSO 

5.0  40  (32dB) 

0 7m 
5.0  _  0.20 


0.55  =  4.4  ( 12.9dB) 

0.55  ,0.07 
ft  n 


0.75  _  2.5  (8dB) 


LOAD  2 


LOAD  1 
LOAD  2 


LOAD  1 
LOAD  2 


Note:  For  steel,  all  numerators  come  from  Fig.  16,  and  all  denominators  come 
from  Figs.  28  and  29  for  Loads  1  and  2,  respectively. 

For  aluminum,  all  numerators  come  from  Fig.  22,  and  all  denominators  from 
Fig.  30  (Load  1)  and  31  (Load  2).  The  numbers  in  parentheses  are  20  log  the 
value  of  the  ratio,  next  to  them. 


It  remains  to  exhibit  the  spectral  values  of  the  normalized  pressure 
| P3/P0 I  amplitudes  at  selected  points  within  the  shell.  These  are  shown  in 
Figures  32  to  35  for  seven  internal  points  located  in  the  upper-half  of  the 
shell's  interior.  The  coordinates  of  these  seven  points  are  (0,0°),  (b,0°), 
(b,90°),  (b,180°),  (b/2,0°),  (b/2,  90°),  and  (b/2,  180°).  The  peaks  observed 
in  the  spectral  band  which  we  have  displayed  (i.e.,  0<kga<10)  correspond  to 

shel 1 -resonances.  Figure  32  corresponds  to  an  unloaded  steel  shell  of 
thickness  b/a=0.80.  Figures  33,  34  and  35  correspond  to  an  unloaded  steel 
shell  of  thickness  b/a=0.85,  0.90,  and  0.95,  respectively.  All  these  plots 
are  computed  from  Equations  (14)  and  (14').  The  only  set  of  resonances 
displayed  in  these  plots  which  admits  an  immediate  simple  interpretation  are 
the  two  seen  in  the  graph  for  the  shell's  center  (0,0°).  These  resonances  can 
be  shown  to  correspond  to  the  first  two  roots  of  Jo'(k3b)=0,  since  it  is  the 
inner  radius  b  that  bounds  the  approximately  rigid  enclosure  around  the 
internal  air-filled  space.  The  first  two  such  roots  are  3.83  and  7.06,  and 
they  have  to  be  divided  by  the  value  of  the  respective  b/a-ratio  in  each  one 
of  the  graphs.  For  example,  in  Figure  32,  this  pair  of  peaks  are  seen  to 
occur  at  k3a=3. 83/0. 80=4. 78  and  7.06/0.80=8.82  just  as  it  is  theoretically 
predicted.  Analogous  situations  are  observed  in  Figures  33  to  35. 

Figures  36  to  39  are  the  counterparts  of  the  previous  four  figures,  but 
now  the  insonified  steel  shell  is  externally  excited  by  surface  Load  #1.  They 
are  computed  at  the  same  seven  internal  points  within  the  band  0<kga<10,  by 
means  of  Equations  (30),  (26),  and  (36),  and  they  are  plotted  versus  k3a.  The 
persistent  pair  of  resonances  observable  at  the  plot  for  the  shell's  center 
(0,  0°)  are  slightly  shifted  downwards  from  the  corresponding  values  they  had 
in  Figures  32  to  35  due  to  the  fact  that  now  the  shell  is  excited  by  Load 
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#1.  The  shift  is  small  since  Load  #1  just  affects  a  couple  of  the  shell 
modes.  For  any  of  the  other  six  considered  points,  the  resulting  peaks 
correspond  to  numerically  computed  shell  resonances,  which  are  not  easily 
obtainable  or  identifiable  by  analytical  methods.  Each  one  of  these  four 
figures  corresponds  to  each  one  of  the  four  shell  thicknesses  used. 

Figures  40-43  are  the  counterparts  of  Figures  36  to  39,  but  now  for  the 
steel  shell  externally  excited  with  Load  #2.  The  computations  of  the  spectra 
are  carried  on  at  the  same  seven  internal  points  by  means  of  Equations  (30), 

(26),  and  (38),  and  they  are  plotted  versus  k3a.  The  only  set  of  resonances 
admitting  quick  and  easy  interpretations  is  the  strong  pair  in  the  plot  for 
the  shell's  center  (0,  0°).  These  resonances  are  seen  to  shift  downwards  from 
the  values  they  had  in  Figures  32  to  35,  since  now  the  shell  is  excited  by  Load 
#2.  These  spectral  plots  show  the  locations  of  the  resonance  peaks  in  the 
k3a-domain  for  all  the  considered  shell -thicknesses.  All  the  peaks  correspond 
to  numerically  computed  modal  shell  resonances,  since  now  they  are  not  roots 
of  easily  tractable  analytic  expressions.  All  these  modal  fundamentals  and 
overtones  quantitatively  display  the  way  resonances  change  from  one  internal 
point  to  another.  This  completes  the  description  of  the  figures. 


NSWC  TR  85-42 


SECTION  6 
CONCLUSIONS 


We  have  developed  the  exact  formulation  required  to  analytically  and 
computationally  model  the  transmission  of  plane  acoustic  waves  into  submerged, 
elastic,  air-filled,  cylindrical  shells  subjected  to  arbitrary  loads  on  their 
outer  surface.  We  have  generated  the  codes  required  to  calculate  these 
interior  fields  in  a  variety  of  pertinent  cases.  The  surface  loads  were 
Fourier-decomposed,  and  in  view  of  the  linearity  of  the  system  of  partial 
differential  equations  governing  this  problem,  we  were  able  to  use  the 
superposition  principle  to  separate  the  effect  of  the  incident  wave  from  the 
effect  of  the  surface  loads  on  the  internally  transmitted  pressure  fields. 
Modal  shell  resonances  were  caused  by  both  types  of  excitations  and  we 
presented  a  study  of  their  interplay  as  they  influenced  the  interior  fields. 
The  exact  full  equations  of  (2-0)  elastodynamics  were  used  to  describe  the 
shell  motions,  thus  going  beyond  any  description  that  could  possibly  be 
provided  by  any  (approximate)  shell  theory.  We  used  the  separated 
contributions  o  the  internal  field  caused  by  the  surface  load  and  by  the 
incident  wave,  to  compute  an  advantage  ratio  which  measured  the  gains  that 
resulted  from  sensing  the  internally  transmitted  fields  rather  than  the 
pressure  fields  outside  the  shell  on  its  external  surface.  This  ratio  showed 
gains:  (1)  that  were  always  larger  than  one,  (2)  that  at  the  internal 
caustics,  present  at  any  frequency,  were  many  times  larger  than  one  (see 
Tables  1  and  2)  ,  and  (3)  that  in  the  vicinity  of  certain  prominent  and 
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persistent  shell  resonance-frequencies,  could  exceed  many  times  over,  the 
results  shown  in  Tables  1  or  2.  The  internal  fields  were  displayed  by  means 
of  isobaric  contours  that  clearly  depicted  the  development  of  the  internal 
field  and  of  its  caustics  as  the  k3a  values  were  increased.  Special  codes 
were  developed  to  generate  these  isobaric  contour  plots  for  all  the  pertinent 
cases  (i.e..  Figures  3  to  6)  which  included:  (1)  the  surface  loads  alone 
acting  on  the  shell,  (2)  the  incident  wave  acting  alone  on  the  shell  and, 

(3)  the  surface  loads  and  the  incident  wave  both  acting  on  the  shell.  Two 
specific  types  of  surface  loads  were  considered  which  were  harmonic  in  time 
and  symmetric  about  the  direction  of  the  incident  wave.  Load  #1  acted  on  the 
entire  outer  surface  of  the  shell  and  had  central  symmetry.  Load  #2  acted 
only  on  a  portion  of  the  shell's  surface  and  had  no  such  symmetry.  In  all 
cases  and  for  all  loads  we  generated  a  succession  of  plots  for  four  increasing 
k3a  values  in  the  band  0  <  k3a  <  10.  The  internal  SPL  was  also  studied, 
computed,  and  displayed  at  (several)  points  inside  the  upper  half  of  the 
shell,  as  a  function  of  k3a.  These  spectral  plots  were  produced  for  various 
loaded  and  unloaded,  insonified  or  silent  combination  of  conditions.  At  all 
the  locations,  these  spectral  plots,  which  accounted  for  all  the  modal 
contributions  to  the  internal  field,  quantitati vely  showed  (viz..  Figures  7  to 
9)  the  distorting  effect  the  surface  loads  had  on  the  spectral  characteristics 
of  the  internal  part  of  the  field  caused  by  the  incident  wave.  The  isolated 
effect  of  the  monopole  mode  (n  =  0)  was  studied  separately  (Figures  10  to  15) 
for  the  three  possible  combinations  of  loading  and  insonification 
conditions.  The  internal  SPL  was  plotted  for  that  n  =  0  mode  alone  as  a 
function  of  k3a  at  various  internal  points  (Figures  10  to  12),  or  as  a 
function  of  radial  coordinate  at  fixed  values  of  k3a  (Figures  13  to 
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15).  In  many  cases  we  showed  which  resonance  peaks  were  caused  by  which  type 
of  excitation  mechanism.  In  some  cases  we  demonstrated  that  the  locations  of 
certain  resonance  peaks  or  resonance  nulls  were  governed  by  very  simple 
resonance  conditions  that  one  may  have  expected  to  remain  veiled  behind  the 
rather  formidable  general  and  exact  formulation  presented  here.  This  was 
particularly  striking  for  the  peaks  in  the  total  pressure  spectrum  at  the 
center  (r  =  0)  of  the  shell,  displayed  in  Figure  7,  and  for  the  nulls  in  the 
radial  dependence  of  the  monopole  contribution  shown  in  Figures  13  to  15. 
These  welcomed  results  served  as  natural  check  points.  In  all  the  results 
presented  here,  the  surface  loads  had  amplitudes  no  larger  than  that  of  the 
incident  wave.  This  meant  that  the  external  “signal -to-noise"  ratios 
considered  here  were  never  less  than  unity.  The  preceding  systematic 
investigation  of  the  fields  transmitted  inside  loaded  shells  at  low-to-medium 
frequencies  showed  clearly  the  filtering  and  focusing  effects  of  the  shell  in 
the  frequency  and  the  spatial  domains,  respectively.  Thus,  for  not  too  over¬ 
powering  suface  loads,  the  thin  air-filled  shell  studied  here  acted  as  a  very 
effective  sound  concentrator,  particularly  within  certain  narrow  spectral 
windows  determined  by  the  resonances  of  the  monopole  vibrational  mode. 
Furthermore,  the  resulting  caustics  within  the  displayed  isobaric  countours 
did  not  seem  to  migrate  too  much  from  their  locations  in  the  unloaded 
configuration,  by  the  distorting  action  of  either  type  of  surface  load.  This 
permanence  of  the  caustic  locations  at  a  given  frequency,  under  various  types 
of  external  loads,  is  a  very  encouraging  result. 

The  final  section  studied  the  effects  of  shell  thickness  and  material 
composition  on  sound  transmissibility  into  the  structure.  Four  shell 
thicknesses  spanning  a  wide  range  of  interest  were  used,  as  well  as  a  couple 
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of  pertinent  structural  materials.  For  a  fixed,  relatively  low  frequency  of 
the  incident  wave,  for  which  the  internally  transmitted  pressure-field  already 
showed  sufficient  development,  we  generated  the  contour  plots  of  the  internal 
SPL  for  all  shell  thicknesses  and  materials.  This  was  done  for  the  insonified 
but  otherwise  unloaded  shell  in  Figures  16  and  22.  When  the  shell  was 
externally  excited  in  addition  by  Load  #1,  the  results  were  shown  in  Figures 
17  and  23.  The  case  of  the  submerged  shell  loaded  by  the  unsymmetric  Load  #2 
was  studied  in  Figures  18  and  24.  For  a  higher  frequency  value  of  k3a=7, 
which  coincides  very  closely  with  the  k3a-value  corresponding  to  the  spectral 
monopol e  resonance  of  the  shell  at  its  center,  we  repeated  the  above 
calculations.  The  results  were  shown  in  Figures  19  to  21  and  25  to  27. 

Internal  SPL  for  shells  that  are  excited  by  the  external  loads  but  not  by 
the  incident  wave,  were  generated  via  the  superposition  principle  and 
displayed  in  Figures  28  to  31  for  all  relevant  shell -thicknesses  and 
compositions.  These  figures  and  also  Figures  16  and  22  permit  the  computation 
of  the  advantage  ratio,  PR,  defined  in  Equation  (31),  for  all  shell 
thicknesses  and  loading  conditions  as  shown  in  Table  2.  Spectral  responses  of 
the  internal  field  at  selected  inner  points  and  for  the  four  shell -thicknesses 

in  question,  were  computed  and  displayed  in  the  final  set  of  12  graphs.  The 

first  four  figures  of  that  set  (i.e..  Figures  32  to  35)  are  for  the  unloaded 
shell.  The  next  four  (i.e..  Figures  36  to  39)  are  for  the  shell  excited  by 
Load  #1  and  the  final  four  (i.e..  Figures  40  to  43)  are  for  the  shell  excited 

by  Load  #2.  This  analysis  of  shell  thicknesses  and  composition  effects  on  the 

internally  transmitted  fields  shows  that  under  the  conditions  of  the  study, 
the  earlier  conclusions  continue  to  hold.  Therefore,  in  view  of  the  lens-like 
sound-concentrating  effect  of  the  shell  in  space  and  of  its  filtering  action 
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in  frequency,  it  appears  more  advantageous  to  sense  internal  rather  than 
external  fields  for  all  the  shell  thicknesses  and  compositions  considered 
here.  Ideal  locations  to  perform  the  internal  sensing  have  been  determined  by 
means  of  the  local  maxima  (i.e.,  the  caustics)  that  we  have  extensively 
computed  in  this  idealized  model  to  quantitatively  illustrate  and  demonstrate 
the  concept. 

For  most  cases  where  this  analysis  will  be  of  practical  value,  the 
assumptions  used  here  seem  to  be  all  justified.  The  incident  waveform  usually 
has  a  fixed  frequency  relative  to  the  target  size.  This  frequency  is  usually 
low.  When  the  incident  frequency  is  high,  the  internal  fields  are  very 
complicated  but  then,  other  mechanisms  enter  into  play  and  this  type  of  study 
loses  its  urgency.  The  incident  wave  and  the  surface  excitations  both  induce 
target  eigenfrequencies  whose  effects  become  manifest  in  the  plots  of  the 
internally  transmitted  fields.  These  resonance  features  became  obvious  when 
plots  are  generated  versus  frequency  or  in  the  space  domain  by  means  of  the 
contour  plots.  It  is  well  known  and  amply  documented  in  the  literature  of  the 
r$t19“24  anai0g0us  resonance  features  manifest  themselves  in  the 
scattered  fields  present  in  the  outer  medium.  The  external  "signal -to-noise" 
ratios  are  generally  high  and  never  fall  below  unity.  The  shell  vibrations 
are  sufficiently  small  to  be  correctly  modeled  by  the  linear  equations  of  two- 
dimensional  axisymmetric  elastodynamics.  Finally,  although  surface 
excitations  caused  by  flow  are  extremely  complicated,  time-varying,  and  in 
some  instances,  probabilistic  in  their  description,  in  specific  cases  they  can 
be  modeled  by  the  simple  types  of  symmetric  or  asymmetric  loads  displayed  in 
Figures  1  and  2. 
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Portions  of  this  work  were  presented  at  a  Professional  Society  Meeting,' 

and  highlights  of  the  most  salient  features  were  published  in  a  Journal 
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Article.  Conversations  with  the  sponsors  at  ONR  (Code  432)  were 
particularly  enlightening  and  we  acknowledge  them  here. 


FIGURE  1.  THE  CYLINDRICAL  SHELL,  WITH  LOAD  NO.  1  ACTING  ON  ITS  ENTIRE  OUTER 
SURFACE,  AS  IT  INTERACTS  WITH  THE  INCIDENT  SOUND  WAVE  OF 
AMPLITUDE  p0 
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FIGURE  4.  INTERNALLY  TRANSMITTED  ISOBARIC  CONTOUR  PLOTS  OF  AN  AIR-FILLED 

CYLINDRICAL  STEEL  SHELL,  LOADED  WITH  LOAD  NO.  1 ,  WHICH  INTERACTS  WITH 
A  PLANE  INCIDENT  ACOUSTIC  WAVE  AT  FREQUENCIES  SUCH  THAT  k,a  -  1, 3, 6  AND  10 
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FIGURE  6. 


INTERNALLY  TRANSMITTED  ISOBARIC  CONTOUR  PLOTS  FOR  A  SUBMERGED 
AIR-FILLED,  CYLINDRICAL  SHELL,  LOADED  WITH  LOAD  NO.  1  (UPPER  PLOT), 
OR  WITH  LOAD  NO.  2  (LOWER  PLOT),  BOTH  FOR  k3a  -  5.  NO  INCIDENT  WAVE 
INTERACTS  WITH  THE  SHELL  HERE 
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FIGURE  9a.  INTERIOR  PRESSURE  LEVELS  P3/Pq  VERSUS  k~a  AT  SIX  (OF  ELEVEN) 
POINTS  INSIDE  A  SUBMERGED  CYLINDRICAL  SHELL,  LOADED  WITH 
LOAD  NO.  2,  AND  INSONIFIED  BY  AN  INCIDENT  PLANE  WAVE 


r 


NSWC  TR  8542 


» 

t 


P(r)i 


« 


NSWC  TR  8642 


FIGURE  11. 


INTERIOR  SPL,  |  P3/P0  | ,  VERSUS  k3«,  COMPUTED  AT  SIX  RADIAL  LOCATIONS 
INSIDE  A  CYLINDRICAL  SHELL  LOADED  WITH  LOAD  NO.  1,  FOR  ITS 
MONOPOLE  <n  -  0)  MODE 
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FIGURE  15. 


RADIAL  DEPENDENCE  OF  THE  INTERNAL  SPL,  |  P3/P0 1 ,  FOR  THE  MONOPOLE 
(n  -  0)  MODE  OF  A  CYLINDRICAL  8HELL  LOADED  WITH  LOAD  NO.  2,  AT  FOUR 


FREQUENCIES  SUCH  THAT  kgt  -  1, 3,  B  AND  10 


CAUSTICS  ARE  DETERMINED 


FIGURE  16.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSl 
LEVELS  INSIDE  AN  UNLOADED  STEEL  SHELL  IN  WATER  AT  kga  -  5  FOR  FOUR 
SHELL-THICKNESSES  (i.a.,  b/a  -  0.80, 0.85, 0.90  AND  0.95).  THE  LOCATIONS  OF  THE 
CAUSTICS  ARE  DETERMINED 


FIGURE  17.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE 
LEVELS  INSIDE  A  STEEL  SHELL  IN  WATER  LOADED  WITH  LOAD  NO.  1  AT  k3a  -  5,  FOR 
FOUR  SHELL-THICKNESSES  (i.a.,  b/a  -  0.80,  0.86, 0.90  AND  0.96) 


FIGURE  16.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE 
LEVELS  INSIDE  AN  UNLOADED  STEEL  SHELL  IN  WATER  AT  k3a  -  5  FOR  FOUR 
SHELL-THICKNESSES  (i.e.,  b/a  -  0.80, 0.85, 0.90  AND  0.95).  THE  LOCATIONS  OF  THE 
CAUSTICS  ARE  DETERMINED 


-JE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE 
WATER  LOADED  WITH  LOAD  NO.  1  AT  k,a  -  5,  FOR 


FIGURE  18.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE 

LEVELS  INSIDE  A  STEEL  SHELL  IN  WATER  LOADED  WITH  LOAD  NO.  2  AT  k^a  ■  6,  FOR  FOUR 
SHELL-THICKNESSES  (i.a.,  b/a  -  030, 0.86, 0.90  AND  0.96) 


FIGURE  20.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE 
LEVELS  INSIDE  A  STEEL  SHELL  IN  WATER  LOADED  WITH  LOAD  NO.  1 ,  AT  k3a  •  7.0, 
FOR  FOUR  SHELL-THICKNESSES  (i.a.,  b/a  -  0.80, 0.85, 0.90  AND  0.95) 


FIGURE  21.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE 
LEVELS  INSIDE  A  STEEL  SHELL  IN  WATER  LOADED  WITH  LOAD  NO.  2  AT  kga  -  7.0,  FOR 
FOUR  SHELL-THICKNESSES  (i.e..  b/a  =  0.80, 0.85, 0.90  AND  0.95) 


FIGURE  22.  ISOBARIC  CONTOUR  PLO  i*  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE 

LEVELS  INSIDE  AN  UNLOADED  ALUMINUM  SHELL  IN  WATER  AT  k,a  -  5,  FOR  FOUR  SHELL- 
THICKNESSES  {i.e.,  b/a  -  0.80, 0.85, 0.90  AND  0.9"  "HE  LOCATIONS  OF  THE  CAUSTICS 
ARE  DETERMINED 


FIGURE  23.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE  LEVELS 
INSIDE  AN  ALUMINUM  SHELL  IN  WATER  LOADED  WITH  LOAD  NO.  1  AT  k3«  -  5,  FOR  FOUR 
SHE LL'THICKNESSES  (i.e.,  b/a  -  0.80,  OSS,  0.90  AND  0.96) 


FIGURE  26.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE  LEVELS 
INSIDE  AN  ALUMINUM  SHELL  IN  WATER  LOADED  WITH  LOAD  NO.  1  AT  k3a  =  7.0,  FOR  FOUR 
SHELL-THICKNESSES  (i.e.,  b/a  =  0.80, 0.85, 0.90  AND  0.95) 


FIGURE  27.  ISOBARIC  CONTOUR  PLOTS  FOR  THE  INTERNALLY  TRANSMITTED  SOUND  PRESSURE  LEVELS 
INSIDE  AN  ALUMINUM  SHELL  IN  WATER  LOADED  WITH  LOAD  NO.  2  AT  k3a  -  7.0,  FOR  FOUR 
SHELL-THICKNESSES  (i.e.,  b/a  =  0.80, 0.85, 0.90  AND  0.95) 


If 


FIGURE  31.  INTERNALLY  TRANSMITTED  ISOBARIC  CONTOUR  PLOTS  FOR  A  SUBMERGED,  AIR-FILLED 
ALUMINUM  SHELL  EXTERNALLY  EXCITED  BY  LOAD  NO.  2,  AT  kga  -  5.0.  NO  INCIDENT 
WAVE  IS  PRESENT  HERE.  FOUR  SHELL-THICKNESSES  ARE  DISPLAYED  HERE:  b/«  -  0.80, 0.85, 
0.90  AND  0.95.  THIS  IS  THE  EFFECT  ON  THE  SURFACE  LOAD  ALONE  AS  GIVEN  IN  EQUATIONS 
(30)  AND  (26) 
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FIGURE  32. 


INTERIOR  PRESSURE  LEVELS  |  p,/p JaT  SEVEN  POINTS  INSIDE  AN  UNLOADED  STEEL 
SHELL  IN  WATER  INSONIFIED  BY  A  PLANE  WAVE,  VS.  k3a.  THE  RELATIVE  SHELL- 
THICKNESS  IS  b/a  -  0.80 
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FIGURE  35.  INTERIOR  PRESSURE  LEVELS  |  P3/pJaT  SEVEN  POINTS  INSIDE  AN  UNLOADED  STEEL 
SHELL  IN  WATER  INSONIFIED  BY  A  PLANE  WAVE,  VS.  k3a.  THE  RELATIVE  SHELL- 
THICKNESS  IS  b/a  -  0.95 
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APPENDIX  A 

NON-VANISHING  ELEMENTS  OF  THE  MATRICES 

The  30  non-vanishing  elements  of  the  matrices  D  and  A*  are  as  follows 
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